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I. INTRODUCTION 



Determination of interface structure from X-Ray diffraction data, involves extracting the 
relative phases between the structure factors of the different components that make up the 
interface. Atomistic modeling of the so-called Crystal Truncation Rods (CTR's)jl] allowed 
the determination of several interfaces, as well as serving as a probe for the strain-induced 
polarization that emerges in a deposited thin filmQ. Electron density (ED) modeling, 
which does not rely on atom positions, allow in principle to search the parameter space of 
the fitted interface model in a more effective way. Techniques to make the ED model to 
converge faster have appeared recently, and of particular note is COBRA^. COBRA has 
the advantage of convergence within a small number of iterations, and was succesfully used 
to determine the structure of Gadolinium Oxide film on GaAs 41, and more recently the 

n 

problem of stress-induced polarization was also tackled [5fl. 

However, ED methods generally steer clear of anomalous scattering. This is done in 
order to simplify the model - by ignoring anomalous scattering, the electron charge density 
is taken to be a purely real spatial function p( r ), and its Fourier transform obeys /( q ) = 
f(—~q)*- This relation allows us to iteratively correct the Fourier transform (or Complex 
Scattering Factor in COBRA terminology), both for positive and negative values of out-of- 
plane momentum transfer (i.e. +q and -q along each Bragg rod), without having measured 
the negative part of the rod. 

Including anomalous scattering in our analysis has the following advantages: 

1. Greater accuracy when heavy atoms are modeled. For example, the magnitude of the 
imaginary correction to the atomic structure factor of Sr[6!] is larger than 1 electron 
almost across the entire energy range 5-20keV. 

2. The anomalous scattering correction acts as a chemical signature, and allows to dif- 
ferentiate between atoms of similar atomic numbers. 

3. In the method shown here, the changes made to the CSF take into account 4 exper- 
imental points per q, instead of 2 as in regular COBRA. Thus reducing the effect of 
spurious measurements. 

We now describe a method by which measurements in two different energies, allows solving 
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for f(— q), without constraining the charge density to be real everywhere. An application 
of this method will appear in a future publication. 

II. METHOD 

Assume that the scattering intensity was measured for two energies, denoted a and b. 
Also, assume only one chemical element of the scattering atoms experiences a subtantial 
change of atomic scattering factor between the two energies. In this situation there is a 
definite relation between the CSF's measured at both energies: 

hCq) — Ja{~q) — AF • (feiCq)) 

Where AF is the difference in the anomalous scattering correction between the two ener- 
gies, and fei{~q) is the geometrical structure factor of the chemical element that contributes 
the most to AF. (i.e. it is the Fourier transform of the spatial distribution of this element). 
~q denotes location along a rod, we can write the same equation for the negative side of the 
rod: 

/rt-/ a (^)=AF.(SH)) 

As noted, f e i{q) is the transform of a spatial density, which is a real function. Therefore 
fe\(~q) = fe\(—~q)*, and we can write: 

(f b (~q) ~ U t ))*exp(2*arg(AF)) = £(=}) - £(=}) 

For a small violation of Friedel's law we can say that |/ (— q)\ = \f a (~q)\, an d now we 
can solve for %(-q), f a (-q) 

K(=?) + (fb(t) - / a (t))*exp(2zarg(AF))| = \f b tf)\ 

So after a new value is arrived at for fb(~q), /a(~<f)by using the usual COBRA iterative 
scheme, we can solve for the negative CSF's. Usually the equation will yield two solutions, 
and using the same arguments as in regular COBRA, we choose pairs of solutions in which 
the correction to the CSF changes slower between consecutive wavevectors. 
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